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ABSTRACT 

We have performed deep imaging surveys for Lyman a emitters (LAEs) at redshift ^7.3 in two blank 
fields, the Subaru Deep Field (SDF) and the Subaru/XMM-Newton Deep survey Field (SXDF), using the 
Subaru/Suprime-Cam equipped with new red-sensitive CCDs and a new narrow-band filter, NB1006 (A c = 
10052A, FWHM AA = 214A). We identified four objects as LAE candidates that exhibit luminosity ex- 
cess in NB1006. By carrying out deep follow-up spectroscopy for three of them using Subaru/FOCAS and 
Keck/DEIMOS, a definitively asymmetric emission line is detected for one of them, SXDF-NB 1006-2. As- 
suming this line is Lya, this object is a LAE at z = 7.215 which has luminosity of 1.2^g x 10 43 erg s" 1 and a 
weighted skewness S w = 4.90 ±0.86. Another object, SDF-NB 1006-2, shows variable photometry and is thus 
probably a quasar (QSO) or an active galactic nucleus (AGN). It shows an asymmetric emission line at 10076A 
which maybe due to either Lya at z = 7.288 or [O II] at z = 1.703. The third object, SDF-NB1006-1, is likely 
a galaxy with temporal luminosity enhancement associated with a supernova explosion, as the brightness of 
this object varies between the observed epochs. Its spectrum does not show any emission lines. The inferred 
decrease in the number density of LAEs toward higher redshift is nf^ 7 a 3 / n Z j~ 5 a 7 = 0.05^^1 from z = 5.7 to 7.3 

down to L Lya = 1.0 x 10 43 erg s" 1 . The present result is consistent with the interpretation in previous studies 
that the neutral hydrogen fraction is rapidly increasing from z = 5.7 to 7.3. 

Subject headings: cosmology: observations — early universe — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Cosmic reionization is one of the most important events in 
the early history of the universe. It is the drastic phase transi- 
tion of the intergalactic hydrogen medium (IGM) from neutral 
to ionized; therefore it is closely related to the formation and 
evolution of the first stars and galaxies. How and when this 
cosmic reionization process took place remains an issue of 
great interest in the studies of the early history of the universe 
and many observational studies have been made to probe 
early galaxies during the last decade (lye 2011). The latest 
measurements of the polarization of the cosmic microwave 
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background (CMB) by the Wilkinson Microwave Anisotropy 
Probe (WMAP) constrained the optical depth of electron scat- 
tering during reionization, suggesting that the probable red- 
shift of reionization, if it took place instantaneously, was 
Z = 1 1 .0 ± 1 .4 (Dunkley et al. 2009). Measurements of Gunn- 
Peterson (GP) troughs (Gunn & Peterson 1965) on spectra of 
high-z quasars sampled by Sloan Digital Sky Survey (SDSS) 
imply that reionization had been completed by z ~ 6, when 
the fraction of IGM neutral hydrogen is estimated to be very 
small, x 1 ^ 6 2 ~ 0.01 - 0.04 (Fan et al. 2006), where x HI is 
defined as xhi = tiHi/riH with rim and n# being the number 
density of neutral hydrogen and all hydrogen in the universe, 
respectively. 

Another promising technique to probe the reionization pro- 
cess is to use Lyman a emitters (LAEs). Narrow-band (NB) 
(bandwidths of ~ 100 -200 A) imaging can efficiently isolate 
LAEs, which prominently emit Lya from their interstellar gas 
ionized by massive stars. Comparing Lya luminosity func- 
tions (LFs) of LAEs before and after cosmic reionization pro- 
vides an approach to determine the epoch of end of reion- 
ization (Haiman & Spaans 1999; Malhotra & Rhoads 2004), 
because Lyman a photons are sensitive to neutral IGM abun- 
dances. In the reionization era, Lyman a photons would be 
scattered at higher redshifts due to the damping wing from the 
surrounding neutral IGM. Therefore, the observed abundance 
of LAEs at this epoch should vary with the neutral fraction of 
IGM hydrogen. Using LAEs to probe reionization is sensitive 
to neutral IGM fractions xhi > 10~ 3 . This is an advantage over 
the GP test which is not sensitive at xhi > 10~ 3 . LAE imaging 
surveys are a good observational method to probe the epoch of 
the cosmic reionization, although it is difficult to distinguish 
the effects of reionization from those of galaxy evolution. 

Due to their very low number density, wide field deep sur- 
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vey observations are essential to construct a statistical sam- 
ple of high-redshift LAEs. Wide field coverage also allows 
us to overcome the uncertainties of cosmic variance, as well 
as potentially large cosmic variance. The Subaru Prime Fo- 
cus Camera (Miyazaki et al. 2002, Suprime-Cam) has ten 
MIT/LL 2048 pix x4096 pix CCDs, located at the prime fo- 
cus of Subaru Telescope, covering a 34' x 27' field-of-view 
with a pixel scale of 0".202. The CCDs of Suprime-Cam 
were replaced to new fully depleted CCDs made by Hama- 
matsu Photonics in 2008 (Kamata et al. 2008) (see the next 
section). Because the practical ranges of redshifts observable 
for LAEs are limited by OH airglow lines, some dedicated NB 
filters, corresponding to z = 3.1,4.8,5.7 and 6.6 were manu- 
factured (e.g., Kashikawa et al. 2006, 201 1; Ouchi et al. 2008, 
2010; Shimasaku et al. 2006; Taniguchi et al. 2005). Sur- 
veys performed using these NB filters revealed that the Lya 
LF exhibits almost no evolution from z = 3.1 to 5.7, whereas 
the Lya LF at z = 6.6 shows an apparent deficit, at least at 
its bright end, compared to that observed at z = 5.7, while the 
rest-UV LFs at z = 6.6 and 5.7 are comparable. This result 
strongly supports the suggestion that the decline in the LAE 
luminosity function at z > 6 was caused by the attenuation 
of the IGM rather than by the galaxy evolution. The derived 
L* at z = 6.6 is 30% -60% lower than that at z = 5.7, indi- 
cating the luminosity attenuation toward higher redshift. On 
the other hand, Hu et al. (2010) reported that L* is unchanged 
from z = 5.7 to 6.5 based on their spectroscopically sample. 
Then, when and how the end of cosmic reionization was com- 
pleted is still under debate. Similarly, discovery of an LAE at 
z ~ 7.0, IOK-1 at z = 6.96, was made using the NB973 filter 
(A c = 9755A, AA = 200A) imaging (lye et al. 2006; Ota et al. 
2008), revealing that the Lya luminosity density gradually de- 
creases from z = 5.7 even to 7.0. In fact, the number density 
ofz = 7LAEs is only 18-36% of the density at z = 6.6. These 
results indicate a substantial transition in the IGM ionizing 
state during these early epochs. We might be observing the 
end of the cosmic reionization process. 

To trace the neutral fraction as a function of redshift, it 
is necessary to survey for LAEs at higher-z. LAEs beyond 
Z = 7 have been surveyed by NIR instruments covering wide 
field-of-view and equipped with NB filters including the 4 m 
Mayall Telescope/NEWFIRM (Hibon et al. 2010; Tilvi et al. 
2010; Krug et al. 2012). They reported that the Lya LF has 
not evolved beyond z = 6. These photometric LAE candi- 
dates have not yet been spectroscopically confirmed and pos- 
sible presence of remaining contaminants cannot be ruled out. 
Therefore, photometric constraining of LFs at z > 7 is not yet 
robust, and spectroscopic confirmation of the Lya line is re- 
quired. 

Another effective method is to measure the fraction of Lya 
emitting galaxies among LBGs in order to understand the 
reionization process. Sizable LBG samples at z ~ 7 found by 
the HSTAVFC3 have been confirmed spectroscopically with 
the Keck/DEIMOS and the VLT/FORS2 (Fontana et al. 2010; 
Ono et al. 2012; Pentericci et al. 2011; Schenker et al. 2012; 
Vanzella et al. 201 1). Their results suggest that the Lya frac- 
tion at z ~ 7 could decline from z < 6 (Stark et al. 2010, 201 1) 
and it might be caused by the increased neutral fraction at 
z ~ 7. Hence, the current suggested redshift range of the end 
of the cosmic reionization from <~ 6 to ~ 11 indicates that the 
details of the cosmic reionization process remain obscure. 

In this paper, we present the results of photometric and 
spectroscopic surveys for z = 7.3 LAEs and use this sample 



to constrain the history of cosmic reionization. The remain- 
der of this paper is organized as follows. First we describe 
the imaging survey for z = 7.3 LAEs in two blank fields in 
§2. We also explain the selection criteria of z = 7.3 LAE can- 
didates based on narrow-band and broad-band photometry in 
§3. Next, we present deep spectroscopic observations of our 
photometrically-selected LAE candidates in §4. We compare 
the estimated Lya luminosity and SFR of z = 7.3 LAEs with 
those of z = 6.6 and 5.7 LAEs and discuss constraints on cos- 
mic reionization in §5. In the last section, we summarize our 
results. 

Throughout the present paper, we adopt concordance cos- 
mology with (fi m , fi A ,/t) = (0.3, 0.7,0.7), (Spergel et al. 2007) 
and AB magnitudes with 2" diameter apertures for the 
NB1006 surveys unless otherwise specified. 

2. IMAGING OBSERVATIONS 

The sensitivity at l/j,m of the previous Suprime-Cam was 
only about 20%. The CCDs of Suprime-Cam were upgraded 
to new red-sensitive ones in 2008 (Kamata et al. 2008). As 
shown in Figure 1 , these new CCDs have a quantum efficiency 
of 40% at twice that of the previous CCDs and enable an 
effective survey for LAEs beyond z = 7. Actually, the survey 
using the new CCDs for LAEs at z = 7 have obtained a deeper 
limiting magnitude (Ota et al. 2010a) than that attained by the 
previous CCDs with same integration time. 

We built a new NB filter, NB1006, to carry out a LAE 
survey at z = 7.3 with the new red-sensitive Suprime-Cam. 
The transmission band of NB1006 filter is centered at A c = 
10052A, and has a FWHM AA = 214A to cover 9943- 
10142A, a modest, but not perfectly clean, OH band gap. 
LAEs in the redshift range z = 7.261 ± 0.082 can be detected 
with this filter (see Figure 1). 

We observed two sky regions, the Subaru Deep Field 
(SDF) and the Subaru XMM/Newton Deep survey Field 
(SXDF) using the new red-sensitive Suprime-Cam and the 
new filter, NB1006. The SDF is centered at 13 A 24 m 21. ? 4, 
+27^29'23"(72000) (Kashikawa et al. 2004). The field has 
the deepest imaging data from those having a single pointing 
of Suprime-Cam. Deep broadband B,V,R,i',z', and narrow- 
band NB816 (A c = 8160A, AX FWHM = 120A) and NB921 
(A c = 9196A, AX FWH m = 132A) filter images were taken 
(Kashikawa et al. 2004) and are made available for pub- 
lic access. The limiting magnitudes in a 2" aperture at 3a 
are B = 28.45, V = 27.74, R = 27.80, i' = 27 .43, z' = 26.62, 
A/B816 = 26.63, NB921 = 26.54, respectively. All of the im- 
ages were convolved to a common seeing size of 0".98. These 
broad-band and narrow-band images can be downloaded from 
the official website of the SDF project K 

The SXDF is centered at 02 /l 18 m 00. ? 0, -05 d 00'00"(./2000) 
(Furusawa et al. 2008; Ouchi et al. 2008, 2010), and cov- 
ers ~ 1.3 square degree region of sky by five pointings of 
Suprime-Cam (SXDF-Center, North, South, East, West). We 
chose a field to cover the southern half of SXDF-Center and 
the northern half of SXDF-South, where Spitzer SpUDS k , 
SEDS 1 , and HST CANDELS (Grogin et al. 2011; Koeke- 
moer et al. 201 1) data are also available, (see Figure. 2 in Ota 
et al. (2010a)). The limiting magnitudes in the 2" aperture 
at 3cr are B = 28.09, V = 27.78, R = 27.57, i' = 27.62, z' = 

J http://soaps.nao.ac.jp/ 

k http://irsa.ipac.caltech.edu/data/SPITZER/docs/spitzermission/ 
observingprograms/legacy/ 

1 http://www.cfa.harvard.edu/SEDS/ 
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FIG. 1. — Filter transmission curves of the Suprime-Cam broad-band and 
narrow-band filters and the new NB1006 filter are shown with corresponding 
labels. The spiky profile at the bottom represents the OH airglow lines. The 
OH airglow lines are not strong, though not completely absent, at wavelength 
in the NB1006 filter. The long-dashed curve at the top shows the atmospheric 
transmission. The short-dashed and thick solid curves show the quantum 
efficiency of the previous MIT/LL CCDs and the new fully depleted CCDs, 
respectively. 

26.57, 7YB816 = 26.55, NB921 = 26.40 at the SXDF-Center 
and B = 28.33, V = 27.75, R = 27.67, i' = 27.47, z' = 26.39, 
7YB816 = 26.65, NB921 = 26.40 for SXDF-South. For SXDF, 
we matched the PSF of NB 1006 images to the maximum PSF 
size of 0".82 as of the SXDF-Center and South. The limit- 
ing magnitudes of available near infrared bands (/, H, and K) 
for SDF and SXDF are much shallower than those for opti- 
cal bands and we did not use data from J,H,K images for the 
photometric selection of candidates. The transmission curves 
of the optical filters used in the observation of SDF and SXDF 
are shown in Figure 1. Spectroscopic standard stars GD153 
and GD7 1 for SDF and SXDF were imaged during the obser- 
vation to calibrate the photometric zeropoint. 

The SDF and SXDF were observed using the NB1006 filter 
in February - April 2009 and October 2009 in 22.3 and 14.8 
hour integrations, respectively. The observational log is sum- 
marized in Table 1 . We removed some images with poor qual- 
ity from final analysis and the actual integration times used 
in the present paper are 17.2 and 14.3 hours for the SDF and 
SXDF, respectively. After masking out those regions that have 
low S/N ratios, the final effective areas of the SDF and SXDF 
images used were 855 and 863 arcmin 2 , corresponding to the 
survey volumes of (2.96, 2.98) xlO 5 Mpc 3 for the SDF and 
SXDF, respectively. 

3. SELECTION OF PHOTOMETRIC CANDIDATES 

3.1. Photometry 

Source detection and photometry were carried out us- 
ing SExtractor version 2.8.6 (Bertin & Arnouts 1996). We 
searched for objects having an area larger than 5 contiguous 
pixels with a flux greater than 2er excess over the sky sur- 
face brightness. Objects were chosen from the NB 1006 image 
and their photometry was performed on the images in other 
wavebands at the same positions. The double-imaging mode 
of SExtractor was used for this task. The 2"-diameter aper- 
ture magnitudes of detected objects were measured using the 
MAG_APER parameter and the total magnitudes were mea- 
sured using MAG_AUTO. When the LAE candidates were 
selected, the aperture magnitude was used as the magnitude 
of object. Astrometric calibration was performed for images 
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FIG. 2. — The detection completeness of our SDF and SXDF NB1006 im- 
ages per 0.5 mag bin. 

of NB 1006 and the astrometric errors are ~ 0".2. 

Bright stellar halos, saturated CCD blooming, and pixels 
of abnormally high or low flux count spikes were masked 
out in the SDF and SXDF images in all wavebands. The 
low S/N regions at the edge of the image were also masked 
out. The final object catalog constructed from the NB1006 
images contains 31, 168 objects for the SDF and 27,565 ob- 
jects for the SXDF down to NB 1006< 24.83(5ct) for SDF and 
NB1006< 24.57(5cr) for SXDF. 

3.2. Detection Completeness 

Photometric detection completeness in photometry is de- 
fined as the ratio of the number of objects detected by the 
photometry to that of the objects actually present in the uni- 
verse. In order to evaluate the detection completeness of our 
detections down to the limiting magnitude of our NB 1006 im- 
ages, the starlist and mkobjects tasks in IRAF were used. We 
created a sample 20,000 artificial stellar objects having the 
same PSF as on the NB1006 images with a random spatial 
distribution over blank regions of the images and with magni- 
tudes ranging from 21 to 27 mag. The SExtractor was run for 
source detection on these mockup image data in exactly the 
same way as we did in our actual photometry. We used the 
ASSOCiation mode of the SExtractor to avoid sky noise not 
located on the artificial stellar objects. We repeated this pro- 
cedure five times and averaged the obtained detection com- 
pleteness. The results are shown in Figure 2. It is likely that 
the detection completeness exceeds 90% at NB1006 < 23.7, 
and more than 50% even at the 5 sigma limiting magnitudes 
of 7YB1006 = 24.83,24.57, for SDF and SXDF, respectively. 
This detection completeness factors were used to correct the 
number and luminosity densities of z = 7.3 LAEs discussed in 
§5. 

3.3. Selection of z = 7.3 LAE candidates 

We select LAE candidates at z = 7.3 as those objects show- 
ing significant flux in NB1006 but not detected in any of the 
bluer bands. The (z' - 7YB1 006, 7YB 1006) color magnitude dis- 
tribution of all the detected objects in the SDF and SXDF 
are shown in Figures 3 and 4, respectively. The average z' - 
NB1006 color excess in the magnitude range 22 < 7YB1006 < 
24 was about +0.2. The dashed curves in these figures show 
3<t boundaries of z' -NB1006 color distribution defined by the 

equation ±3a BB -m = -2.5 log I0 (l =p (fi aNB + /Lbb) ' 5 I ' fm)- 

As the diagram indicates, if z 1 -NB 1006 > 2.6 is used as the 
color selection criterion, then z = 7.3 LAEs can be identified. 
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TABLE 1 

Summary of the SDF and SXDF Observations 



Field 


[s] 


the number of images 


t b 

l exp 

[s] 


Date of Observations 


Seeing 
[arcsec] 


SDF 


500-1800 


19 


22094.5 


2009 February 25 


0.8-1.1 


SDF 


1200-1800 


15 


22800 


2009 February 26 


0.8-0.9 


SDF 


730-1200 


7 


7931.2 


2009 March 24 


0.8-1.2 


SDF 


1200 


2 


2400 


2009 April 1 


1.4 


SDF 


900-1200 


22 


25200.8 


2009 April 2 


0.8-1.2 


SXDF 


1200 


23 


27600 


2009 October 15 


0.4 


SXDF 


1200 


21 


25200 


2009 October 16 


0.4 



a The integration time of one image. 

b The total integration time of the SDF is 80426.5 s and the total integration time of the SXDF 
is 52800 s. 
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FIG. 3. — z! -NB1006 color as a function of NB1006 magnitude of the 
objects detected in the SDF (shown by dots). The vertical line is the limiting 
magnitude of the SDF, 24.83 mag. The diagonal dashed line is the 3tr limit 
of the z! magnitudes in 2" aperture. The points on this diagonal dashed line 
are the colors which are converted to the 3cr limits of z! magnitudes in 2" 
aperture. The curve shows the 3cr error track of z! -/VB1006 color. The 
squares are the z. = 7.3 LAE candidates, SDF-NB1006-1 and SDF-NB1006- 
2. 

It should be noted that the NB1006 band pass is located at 
the red edge of the CCD sensitive wavelengths; therefore, the 
UV continuum redder than Lya emission is hardly detected 
even in the z'-band. Hence, we applied the following selection 
criteria to the detected objects in the same way as Ota et al. 
(2008): 

1. yVB1006>5o-: 
NB973 7 NB921,NB816,i',R,V,B < 3a,z' -NB10Q6 > 
2.6 

2. /YB1006>5cr: 
NB973,NB921,z' : NB816,i',R,V,B < 3a 

where all magnitudes were measured as 2"-diameter aperture 
magnitudes, MAG_APER. 

Color criterion, z' -NB 1006 > 2.6, was determined by the 
model spectrum of a galaxy at z = 7.3. We created a spectral 
energy distribution (SED) of a starburst galaxy using stellar 
population synthesis model, GALAXEV (Bruzual & Char- 
lot 2003), with a metallicity of Z = Z Q = 0.02, an age of 
t = 1 Gyr, a Salpeter initial mass function with lower and 
upper mass cutoffs of nti = 0.1M Q and ntu = 100M©, and 
an exponentially decaying star formation history for r = 1 
Gyr. These parameters were selected to be the same as those 
used in Taniguchi et al. (2005); Ota et al. (2008), for consis- 
tency. The absorption by IGM was applied to the generated 
spectrum using the model of Madau (1995). The flux of a 
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FIG. 4. — Same as Figure 3, but the filled squares and dots are the candidates 
and objects detected in the SXDF, respectively. The objects in two fields in 
SXDF, SXDF-C and SXDF-S which have different z' limiting magnitudes 
are plotted in this same diagram. The fainter magnitudes of objects than z' 
limiting magnitude of each fields are converted to the shallower 3<r limits 
(SXDF-S). 



Lya emission line with EWq = OA was added to the SED at 
(1 + 7.26)1216A. In other words, we create a SED of LBG 
with no Lya emission line. Hence, the presently adopted color 
criterion, z' -NB1006 > 2.6, corresponds to the rest-frame 
equivalent width of Lya of EWo > 0. Moreover, the photo- 
metric error may cause interlopers at lower-z to happen to sat- 
isfy our color selection criteria. We carried out the following 
test in the same way as did in Bouwens et al. (201 lb); Wilkins 
et al. (2011). First, we choose brighter sources in magni- 
tude range of 22.5 > NB1006 > 23.5 from our two fields and 
then dim these bright sources to match to the same magni- 
tudes as of our LAE candidates by scaling the optical BB and 
NB fluxes. We then distribute these objects over the real im- 
ages. The number of these artificial objects is the same as 
the number of observed objects in the magnitude range from 
7VB1006 < 23.5 down to each NB1006 limiting magnitude. 
We applied exact the same detection and selection criteria as 
the candidate selection. As a result, we found no objects hap- 
pen to come into our criteria, suggesting that our selection 
criteria can robustly avoid any nearby red and line emitting 
interlopers. 

As a consequence of selection, we found no candidates 
meeting only criterion (1) and two objects satisfying both cri- 
teria (1) and (2) in each field. The images and photometric 
properties of these LAE candidates are listed in Figure 5, and 
Table 2. We detected no flux even in the stacked optical and 
NIR images of all the candidates, as shown in the bottoms of 
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Figure 5. 

4. SPECTROSCOPY 

We carried out optical deep spectroscopy for three of our 
four photometric LAE candidates. We used the Faint Object 
Camera and Spectrograph (FOCAS) in Multi-Object Spec- 
troscopy (MOS) mode mounted on the Subaru telescope on 
2010 March 18, 19, and 21 to follow the two candidates in the 
SDF, SDF-NB 1006-1 and SDF-NB 1006-2. The observations 
were made with a 555 mm" 1 grating and an 058 order-cut 
filter, giving spectral coverage of 7500- 10450A with a dis- 
persion of 0.74Apix _1 . The 0".8-wide slit used gave a spec- 
troscopic resolution of 5.7A (R^ 1500) which is sufficient 
to distinguish [O II] doublet lines from possible contaminant 
galaxies at z=1.7. However, it is likely that [O II] doublet 
lines are not resolved in this resolution if the emission lines 
are too faint or are affected by OH sky lines. The observing 
nights were photometric, with good seeing of ~ 0".6- 1".0. 
Because no emission line feature was visible for the brighter 
candidate, SDF-NB 1006-1, in the entire transmission range 
of the NB1006 filter in 5 hours integration, this object was 
abandoned. We switched on to the second candidate, SDF- 
NB 1006-2, and secured spectra for a total integration time of 
11.5 hours. 

In addition to the observations made in March, we observed 
SDF-NB 1006-2 for a half night on 2010 May 5. The see- 
ing 0".8 - 1".0 was somewhat poor but the weather condition 
was photometric. The setup of the instrument was exactly the 
same as that in March. The integration time in May was 4.5 
hours, and total integration time for the SDF-NB 1006-2 spec- 
troscopy is 16 hours. 

Each of the 30 minute exposures was taken by dithering the 
telescope pointing along the slit by ±1".0 under good seeing 
conditions and ±1".5 under moderate seeing conditions to en- 
able more accurate sky subtraction. The standard star Feige 
34 and Hz 44 were taken at the beginning and end of each 
observed night except on May 5 when only Feige 34 was ob- 
served (Massey & Gronwall 1990). Feige 34 and Hz 44 were 
used for flux calibration of May and March data, respectively. 
The object and standard star data were reduced in a standard 
manner with IRAF. Also, we evaluated the slit loss for the Lya 
flux measurements calculated from the spectra. The targets 
were assumed to be almost perfectly centered in the slit and 
point sources. The slit loss was calculated by putting a Gaus- 
sian profile with FWHM derived from the average seeing size 
in the FOCAS spectroscopic observation on the 0".8-wide 
slit. The slit loss of the FOCAS spectroscopy was evaluated to 
be ^23%. The resultant slit loss was corrected for the spec- 
troscopically measured Lya for the object. The stacked one- 
and two- dimensional spectra of SDF-NB 1006-2 are shown in 
Figure 6. A prominent emission line is seen at 10075.8A. 

We also used DEep ImagingMulti-Object Spectrograph 
(DEIMOS) in MOS mode at the Nasmyth focus of the Keck 
II telescope to take spectroscopy for the candidate, SXDF- 
NB 1006-2 as described in (Ota et al. in prep). The seeing 
was - 0".6- 1".0. We used the GG495 filter and the 830 
lines mm" 1 grating, which was tilted to place a central wave- 
length of 9000A on the detectors. This configuration pro- 
vided spectral coverage between 7000 A and 10400A. The 
spatial pixel scale was 0".l 185 pix" 1 , and the spectral disper- 
sion was 0.47Apix _1 . The slit widths were 1".0, which gives 
the spectral resolution of 3. 3 A. The total integration time of 



DEIMOS is 3. 5 hours. 

The reduction of the spectra for SXDF-NB 1006-2 was per- 
formed using spec2d IDL pipeline developed by the DEEP2 
Redshift Survey Team (Davis et al. 2003). The final one- 
dimensional data output from the spec2d, however, show ap- 
parent over-subtraction of the foreground OH sky lines. In or- 
der to solve this difficulty, we modified the code in the process 
of calibrating the wavelength using arc lamp and subtracted 
OH sky lines carefully in a manner similar to that employed 
in the FOCAS data reduction with IRAF. 



Although we also observed the standard star Feige 110 
for flux calibration with DEIMOS, these data were not re- 
duced due to a computational problem. Hence, the flux of 
the SXDF-NB 1006-2 was calibrated using the spectra of five 
late-type stars observed simultaneously in the same slit mask 
mounted on DEIMOS. We compared their spectroscopically 
measured fluxes with those expected from their photometry 
convolved with the filter (NB1006) transmission and derived 
the zero-point for the spectroscopic flux calibration. The sta- 
tistical fluctuation in the above flux calibration with five late- 
type stars was used as the estimated flux error of the SXDF- 
NB1006-2. The slit loss has already been corrected in this 
calibration. The stacked one- and two- dimensional spectrum 
of the SXDF-NB1006-2 are shown in Figure 7. 

We measured the weighted skewness of the detected emis- 
sion lines, which can be used as a criterion to identify high-z 
Lya emission lines, (Kashikawa et al. 2006). The weighted 
skewness was calculated in the wavelength range where the 
flux is higher than 10% of the peak value, ranging from 
10068A to 10088A, and from 9980A to 9998A, for SDF- 
NB 1006-2 and SXDF-NB 1006-2, respectively. They were 
S u = 4.59 ± 1.25 for SDF-NB1006-2 and 4.90 ±0.86 for 
SXDF-NB1006-2. These measured skewness values satisfy 
the usually adopted criterion for the presence of asymmetry 
of high-z Lya emission, S u > 3 for typical LAEs at z > 5.7. 
Note however that the red side of emission line in the SDF- 
NB 1006-2 spectrum has an overlap with a weak OH sky line, 
and the derived skewness value might be overestimated. 

We had no chance to observe the remaining fourth can- 
didate in the SXDF, SXDF-NB 1006-1 during our allocated 
time. The spectroscopic properties of the line emitting objects 
are summarized in Table 4. 

5. DISCUSSION 

5.1. Check for the Variability of the Candidates 

Since the NB1006 imaging observations we used for detect- 
ing our z = 7.3 LAE candidates were made a few years after 
the imaging observations for other bands, some of the objects 
which are only visible in NB1006 could be variable objects, 
such as galaxies with supernovae or enhanced active galactic 
nuclei (AGNs), which happened to increase their luminosi- 
ties during our NB1006 imaging observations. Since our SDF 
imaging was made in three epochs, we created three stacked 
NB1006 SDF images, one for February 2009, one for March 
2009, and one for April 2009 to check for any variability. The 
SXDF observations were taken on two successive nights on 
October 15 and October 16. Although the time interval is 
short for SXDF, we also examined possible photometric vari- 
ability. 

Objects were chosen from the stacked NB1006 image in all 
epochs and photometry was performed on the images in the 
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FIG. 5. — Multi-waveband images of the z = 7.3 LAE candidates in the SDF and the SXDF. The postage stamps have an area of 10" X 10". The optical images 
are the archival data, which were observed in the SDF and the SXDF projects. The near-infrared images were observed by the UKIRT/WFCAM (Hayashi et al. 
2007). Second row from the bottom: Centered-staked postage-stamp images for LAE candidates. Stacks have no flux in the optical- and NIR- filters. Bottom 
row: Centered-stacked postage-stamp images using only optical- and NIR- filters. 



TABLE 2 

Photometric Properties of z = 7.3 Lyo Emitter Candidates 



Name 




RA 




DEC 


z' a 


NB1006 A/>£R 


NB1006 AUT o 






(J2000) 




(J2000) 


[mag] 


[mag] 


[mag] 


SDF-NB 1006-1 


13 


: 24 : 39.993 


+27 


: 39 : 40.67 


> 26.62 


23.96 ±0.03 


23.62 ±0.04 


SDF-NB 1006-2 


13 


: 24 : 35.418 


+27 


: 27 : 27.81 


> 26.62 


24.61 ±0.05 


24.21 ±0.05 b 


SXDF-NB1006-1 


02 


: 18 : 50.657 


-05 


: 24 : 40.25 


> 26.39 


24.45 ± 0.05 


23.81 ±0.05 


SXDF-NB 1006-2 


02 


: 18 : 56.536 


-05 


: 19 : 58.92 


> 26.39 


24.56 ±0.06 


24.26 ±0.07 



a The 3a limiting magnitudes of the candidates in z'-band were derived from the 2"(A-aperture photom- 
etry. 

b The NB1006 magnitude of the SDF-NB1006-2 was derived by 3" diameter aperture photometry 
instead of MAG_AUTO to avoid contamination from the flux of a nearby object. 



TABLE 3 

the Summary of the Follow Up Spectroscopy 



Target 


[s] 


the number of image 


h'xp 

[s] 


Date of Observations 


Instrument 


SDF-NB 1006-1 


1800 


10 


18000 


2010 March 18-21 


Subaru/FOCAS 


SDF-NB 1006-2 


1800 


23 


41400 


2010 March 18-21 


Subaru/FOCAS 


SDF-NB 1006-2 


1800 


9 


16200 


2010 May 5 


Subaru/FOCAS 


SXDF-NB 1006-2 


1800 


7 


12600 


2010 October 10 


Keck/DEIMOS 



a The integration time of one image. 

individual epochs at exactly the same positions as the objects 
in the NB1006 filter using the double-imaging mode of SEx- 
tractor. The results of independent brightness measurements 
at different epochs of observation are summarized in Table 5 
and shown as postage stamp images of the individual epochs 
are shown in Figure 8. 

As shown in the table, the magnitudes of candidates SDF- 
NB1006-1, and 2 vary by more than ~ 0.5 mag over two 
months, while the average magnitude inconsistency of other 
objects in the SDF NB1006 images which have the similar 
magnitudes as the LAE candidates is less than 0. 1 ± 0.6 mag. 



In particular, SDF-NB 1006-1, the brightest among ourz = 7.3 
LAE candidates, appears to show a significant luminosity de- 
cline of 1 .8 mag as shown in Figure 8. The brightness change 
corresponds to la p j wt , where a p k t is the la photometric er- 
ror of the object. This amplitude of variability is roughly 
consistent with the J-band light-curve of a type la supernova, 
depending on the phase of its explosion (Burns et al. 2011). 
We conclude that SDF-NB 1006-1 is a variable object which 
brightened suddenly in February and then gradually faded. It 
may be a galaxy with a supernova, which is consistent with 
the fact that no obvious emission line was confirmed in the 
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TABLE 4 

Spectroscopic Properties of Line Emitting Objects 



Name 


Redshift 


j?Lya.a 

[10~ 17 erg s-'crrr 2 ] 


L Lya 

[10 43 ergs-'] 


SFR Lvob 

[Meyr" 1 ] 


[mag] 


SRF (/Vd 

[Moyr 1 ] 


FWHM e 
[A] 


[A] 


EW f 
[A] 


SDF-NB 1006-2 
SXDF-NB 1006-2 


7.288? g 
7.215 


0.91 ±0.13 
1 Q+2.5 
1,y -09 


0.57±0.15 
1 2 +L5 


5.2 ±1.1 

n +15 
1J --5.3 


-23.79 ±0.04 

-22 4 + °° 


1646±54? h 
4^7+190 


13.4 
11.5 


4.59 ±1.25 
4.90 ±0.86' 


1.99 ±0.37 
15.4^ 



a The observed line flux corresponds to the total amount of the flux within the line profile. A correction was made for slit loss (see §4). 
b The SFRs were derived from L LyOL : SFR Lya = 9. 1 X VS~ t> l}* a M®yr~ i (Kennicutt 1998) 

c The absolute magnitude of the UV continuum were derived from the UV fluxes f JV = /nbioog- f Lya using equation (1) of Kashikawa et al. (201 1). 
d The SFRs were derived from Muv- SFR UV = 1.4 X 10~ n L v (UV)M e yr~' (Kennicutt 1998) 

e The FWHMs (in observed frame) were derived from the profile fitted by the Gaussian function, since the actual spectra are noisy. 

The rest-frame equivalent widths of Lya emission were derived from the combination of f Lya and Muv ■ 
? The NB1006 magnitude of the SDF-NB 1006-2 fluctuated approximately ±0.8 mag, and it might be a QSO/AGN atz = 1.703 or : = 7.288. If SDF-NB 1006-2 is a [Oil] emitter, 
the Lfoii] and SFR IO ii] are 1.77 X 10 41 [erg.? -1 ] and 2.48[Mq}t~'], respectively. 
h This high UV SFR of SDF-NB 1006-2 may originate from the central QSO/AGN engine. 
1 The weighted skewness of the SXDF-NB 1006-2 were derived from the smoothed emission line. 
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FIG. 6. — Top panel: Composite two-dimensional spectrum of SDF- 
NB1006-2. Bottom panel: Combined spectrum of SDF-NB1006-2. The 
dashed line represents the OH sky line. A prominent emission line is seen 
at 10075. 8A corresponding to z, = 7.288 if this line is Lya. The weighted 
skewness of the line is 4.59 ± 1.25. This value indicates that this emission 
line is asymmetric because it is greater than the threshold (Su = 3) used in 
the LAE studies at z = 5.7 and 6.6 (Kashikawa et al. 2006). The S/N ratio of 
the line is ~ 8.3. The dashed line represents OH sky line. 



FOCAS spectroscopy. 

The possible magnitude variation of SDF-NB 1006-2, 
which shows an emission line with an asymmetric profile 
is a bit more enigmatic. The emission line indicates that it 
could be a LAE at z = 7.288 considering its high weighted 
skewness (S u = 4.59 ± 1.25). On the other hand, because the 
NB1006 magnitude of this object changed by ±0.8-0.9 mag 
(3.3-7.9oy,/, ,) over about one month, it could also be an AGN 
or a QSO. This suggestion might be supported by the high UV 
magnitude of SDF-NB 1006-2 derived by subtracting the Lya 
flux from the total flux in the NB1006 magnitude. Consider- 
ing that the line width is narrower than that of typical QSOs, it 
is unlikely to be a QSO. The Lya flux emitted from the broad 
line region (BLR) in AGNs/QSOs varies slightly, while that of 
[O II] irradiated from the narrow line region (NLR) does not 




9900 



9950 



10000 



10050 



obs 



[A] 



FIG. 7.— Top panel: Combined spectrum of the SXDF-NB 1006-2. Bot- 
tom panel: Combined spectrum of the SXDF-NB 1006-2. A line is visually 
identified at 9987. 6A. The S/N ratio of the line is ~ 4.5. Both spectra are 
smoothed by being convolved with the cr = 1 .5 [pixel] Gaussian function. The 
weighted skewness line is 4.90 ±0.86. This value indicates that this emis- 
sion line is asymmetry because it is greater than the threshold (Su = 3) used 
in the LAE studies at z = 5.7 and 6.6 (Kashikawa et al. 2006). The dashed 
line represents OH sky line. 



change (Kaspi et al. 2007). Besides, the change of UV contin- 
uum flux radiated from the BLR is investigated by the struc- 
ture function of the photometric variability of AGNs/QSOs 
(de Vries et al. 2005), and is ~ 0.55 mag for a year. The 
structure function is the ensemble average of QSO variability 
as function of time and has been used as a tool to character- 
ize the QSO variability. Since the UV continuum flux of the 
SDF-NB 1006-1 is only ~ 70% of the NB1006 flux, its large 
variability might not be due to AGNs/QSOs, though the vari- 
ation of the investigated structure function in the literature is 
relatively large. Hence, this object might be an AGN/QSO 
at z = 7.288 if this line is Lya, or at z = 1.703 if that is [O 
II], rather than a LAE, nevertheless the large amplitude of the 
variability cannot be explained adequately. 
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By contrast, the magnitude of the candidates in the SXDF, 
SXDF-NB 1006-1 and 2, did not show any significant vari- 
ation that would lead us to suspect their AGN nature. The 
brightness changes of SXDF-NB 1006-1 and 2 correspond to 
0.2 mag {\-\<J p hot) and 0.42 mag (2.2a p i wt ), respectively. The 
reality of the SXDF-NB 1006-2 at z = 7.215 is also supported 
by this result of checking variability. 

Moreover, we investigate the possibility that the detected 
lines could be Ha, H/3, or [O III] in terms of metallicities es- 
timated from diagnostic line ratios. The same case test for 
z ~ 6.9-7.2 LBGs was also discussed in Ono et al. (2012). 
Any other emission line features were not observed in the 
spectra; therefore the metallicity is only given as an upper 
limit by using the flux limit. If the detected lines were H/3, 
[O III]A5007 line would fall at 10378A and 10288A for 
SDF-NB1006-2 and SXDF-NB1006-2. The upper limits of 
[O ill]A5007/H/3 line ratios are calculated to be < 0.1-0.4, 
corresponding to the metallicities of 12+log(0/H) > 9.0-9.5 
based on the Figure 17 of Nagao et al. (2006). In the case 
of such high metallicities, they should be very massive and 
bright taking account the mass-metallicity relation. However, 
they are not detected in our deep broadband images. If the 
detected lines were [O lll]A5007, the H/3 line would be de- 
tected at 9782A and 9697A for SDF-NB 1006-2 and SXDF- 
NB 1 006-2, and the lower limits of [O III] A5007 /H/3 are given 
as > 1 1 and > 2, respectively. The detected line of the SDF- 
NB 1006-2 is unlikely to be [O III], because no outlier with [O 
m]A5007/H/3 > 10 in Nagao et al. (2006). The deep FOCAS 
spectroscopy for SDF-NB 1006-2 completely rules out the 
possibility of other lines. On the other hand, the lower limits 
for SXDF-NB 1006-2 corresponds to subsolar oxygen abun- 
dances, which may be typical of low-mass galaxies. Hence, 
the non-detection of H/3 cannot strongly rule out the possi- 
bility of the detected line in the SXDF-NB 1006-2 spectrum 
being [O lll]A5007. The possibility of Ha can be ruled out 
because either [O III] or H/3 is expected to be simultaneously 
detected at the observed wavelength for both targets. The pos- 
sibility that the detected line is the unresolved [O II] doublet 
can be ruled out for SXDF-NB 1006-2 due to the spectral res- 
olution of DEIMOS. On the other hand, for SDF-NB 1006-2, 
we cannot deny the possibility of [O II]. 

After all, we could completely rule out the possibility of the 
detected emission lines being other nebular emission lines, i.e. 
Ha, H/3, and [O III]. One possibility that the detected line of 
SXDF-NB 1006-2 is H/3 could not be completely ruled out 
based on the diagnostic line ratios. However, we should note 
that the asymmetric line profile strongly implies the line to be 
Lya emission. The SDF-NB 1006-1 and 2, which have signif- 
icant variability amplitude more than 3a p i wl , are concluded as 
variable objects. We conclude that the SXDF-NB1006-2 is a 
plausible LAE at z = 7.215, but the SDF-NB 1006-2 is likely to 
be a AGN/QSO either at z = 7.288 or at z = 1 .703. Our results 
indicate that checking the time variability is crucial to strictly 
verify the LAE photometric candidates if they have been de- 
tected in only a single band. The checking the variability is 
also attempted in Bouwens et al. (201 la). This check implies 
a possibility that two of the LAE candidates in the SDF could 
be transient objects such as an AGN/QSO or supernova, even 
though an asymmetric line was detected. However, we should 
note that the time span in SXDF was only one day delay which 
was too short to completely trace the typical variability. 

5.2. The Lya Luminosity Function of the z = 7.3 LAEs 




Feb. Mar. Apr. 15, Oct. 16, Oct. 



FIG. 8. — The images of the z = 7.3 LAEs candidates in the divided images 
for each observational epoch. The PSF sizes are matched to 0".98 in SDF 
and 0". 82 in SXDF. 

In Figure 9, we compare the luminosities of our spectro- 
scopically confirmed LAE at z = 7.215 with the previous LF 
results obtained at z < 7.0. The SDF-NB 1006-2 is also in- 
cluded as a preliminary measurement. The number densities 
of both objects are corrected for the detection completeness 
derived in §3.2. The error of the number density of the z = 7.3 
LAEs includes the possible cosmic variance and the Poisso- 
nian error for small number statistics. The values in Table. 

1 and 2 of Gehrels (1986) are used as the upper and lower 
limit of Poissonian errors. The relative cosmic (root) variance 
er v is calculated by cr v = baDM, where b and gum are the bias 
parameter and the (root) variance of dark halo, respectively 
(Somerville et al. 2004). The ctdm at z = 7.3 with our sur- 
vey volume of 5.94 x 10 s Mpc 3 is calculated to be 0.033 by 
the model of Sheth & Tormen (1999). The bias parameter is 
adopted to be b max = 9 derived from the LAE sample at z = 6.6 
detected in the full SXDF, which is ~ 5 times wider than the 
SDF (Ouchi et al. 2010). The resultant cosmic variance a v 
is - 30%. The error of Lya flux of the SXDF-NB 1006-2 in- 
cludes the error in the flux calibration from using the dwarf 
stars (see §4) in addition to the spectral sky noise. The cumu- 
lative Lya LFs derived from the photometric LAE candidates 
at z = 7 and 7.7 (Hu et al. 2010; Hibon et al. 2010, 2012; 
Krug et al. 2012; Tilvi et al. 2010) are also plotted in Fig- 
ure 9. Their photometric Lya LFs at z = 7 and 7.7 in Hibon 
et al. (2010); Ota et al. (2010a) is based on Lya luminosity 
of each LAE candidate directly converted from Lya flux by 
assuming f Lya = f NB . However, Ota et al. (2010b) estimated 
that - 77% and 23% of NB973 total flux of a z = 6.96 LAE 
comes from Lya and UV continuum fluxes, respectively, us- 
ing imaging and spectroscopy data. Hence, in Figure 9 we 
shifted the photometric z = 7 and 7.7 Lya LFs faintward with 
0.77 x L Lya for fair comparison with other spectroscopic Lya 
LFs. As shown in Figure 9, the luminosity of SXDF-NB1006- 

2 is fainter than those of the LFs at z = 5.7 and 6.6, though 
their errors are large. Furthermore, the number density of 
the SXDF-NB 1006-2 is slightly smaller than that of the LAE 
at z = 6.96, IOK-1, which was selected a single pointing of 
NB973 imaging (only SDF) with the Suprime-Cam, while 
their Lya luminosities are similar. What needs to be em- 
phasized is that secure LAE samples can be constructed by 
spectroscopy, line-asymmetry investigation, and checking the 
magnitude variability, and compared appropriately with LFs 
at z < 7.0. While the previous photometric LFs at z = 7 and 
7 .7 do not change significantly from z < 7, the number density 
of the LAE at z = 7.3 shows deficit. 

We also compare the number and SFR densities at z = 7.3 
with those of other LAE surveys at z < 7 down to L Lya = 
1.0 x 10 43 erg s" 1 in Figure 10 in the same way as Ota et al. 
(2008) to investigate the redshift-evolution of these densi- 
ties. The densities of LAEs at z = 5.7 and 6.6 in Kashikawa 
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TABLE 5 

The variability of the z = 7.3 LAE Candidates 



SDF 


Feburuary 
[mag] 


March (Mar. -Feb.) 
[mag] 


April (Apr. -Feb.) 
[mag] 


SDF-NB 1006-1 
SDF-NB 1006-2 


23.67±0.027 
24.53±0.051 


24.31±0.097 (0.64±0.12, 5.3a phot ) 
23.75±0.059 (-0.87±0.11, !S(r pha ) 


25.49±0.23 (1.82±0.26, 1.Qa phat ) 
25.35±0.20 (0.82±0.25, 33a phot ) 


SXDF 


15, October 
[mag] 


16, October (16, Oct. - 15, Oct.) 
[mag] 




SXDF-NB1006-1 

SXDF-NB 1006-2 


24.45±0.055 

24.44±0.054 


24.65±0.12 (0.20±0.18, \\o phot ) 
24.86±0.14 (0.42±0.19, 2.2^,,,) 





NOTE. — All the magnitudes are calculated by 2"-aperture photometry. The aperture magnitudes (2 // ) at each epoch, the 
difference from the magnitude at the first observed epoch, and the variability significance compared with the photometric error 
&piwt are listed. The PSF sizes are matched to '.98 in the SDF and '.82 in the SXDF. 




L(Lya) [erg s ] 

FIG. 9.— Cumulative Lya LFs atz < 7 and the number densities of z = 7.3 LAE, the SXDF-NB 1006-2 (filled square) and the SDF-NB 1006-2 (open square). 
The number densities are corrected using the detection completeness derived in §3.2. The error of the number densities includes cosmic variance and Poissonian 
error. The error of the Lya luminosity includes the spectral sky noise of both objects in the z = 7.3 survey, and the error in the flux calibration of the DEIMOS 
data in addition to the sky noise of SXDF-NB 1006-2 (see §4). The filled circle is the LAE at z = 6.96 identified spectroscopically (Ota et al. 2008). The dashed 
blue and red curves represent the Lya LF fitted by Schechter function of LAEs at z = 5.7 and 6.6 in the SDF, respectively (Kashikawa et al. 201 1). The solid 
blue and red curves are those at z. = 5.7 and 6.6 in the full SXDF, which is ~ 5 wider than the SDF (Ouchi et al. 2010). The dark-green triangles, bold- and 
dashed-curves are photometric z = 7 LAE samples in Ota et al. (2010a), fitted Schechter functions in the D33 and D41 field of the CFHTLS in Hibon et al. (2012) 
derived by using the Suprime-Cam equipped with NB973. The orange filled squares, filled circles, and curve are photometric z = 7.7 LAE samples derived by 
using the Mayall 4m telescope/NEWFIRM camera (Krug et al. 2012; Tilvi et al. 2010; Hibon et al. 2010). The slope parameter a of all the fitted Schechter 
functions are fixed to be -1.5. If SXDF-NB 1006-1 is a real LAE at z = 7.3 with brighter L Lya than one of the SXDF-NB 1006-2, the cumulative number densities 
of SXDF-NB 1006-2 and SDF-NB1006-1 become the small filled square and the open square without L Lya error slightly shifted in the Lya luminosity direction 
for clarification. 



et al. (2006) are updated to those of the ~ 80% spectroscopi- 
cally identified samples (Kashikawa et al. 201 1), and the one 
of LAEs at z = 6.6 selected from the ~ 5 times wider sur- 
vey area (Ouchi et al. 2010) than SDE It is notable that the 
Lya LF of Ouchi et al. (2010) is mainly from the photomet- 
ric LAE sample. The density errors are included as in Fig- 
ure 9. As shown in this figure, compared with number at 
z = 5.7 and 6.6, the LAE number density at z = 7.3 decreases 
by factors of nff a 3 /nf=£ 7 = 0.05!^ from z = 5.7 to 7.3 and 

n % 7 a/ n Ly 6 a 6 = °- 16 -o:i6 from Z = 66 to 7 - 3 ' although the sta- 
tistical uncertainty at z = 7.3 is large. The number and SFR 
density might decrease, nff^/n^ = 0. 72^1 from z = 7.0 



to 7.3. Figure 10 indicates that the number and SFR densities 
monotonically decrease clearly from z = 5.7 to 7.3, while they 
do not change from z = 3 to 6. 

In contrast, the number density of the rest-UV continuum 
of the SXDF-NB 1006-2 at z = 7.3 might be comparable with 
the bright end of the rest-UV continuum cumulative LFs at 
z = 5.7 and 6.6 as shown in Figure 11. These UV continuum 
LFs were constructed from LAE samples in SDF (Kashikawa 
et al. 201 1). The UV magnitude of SXDF-NB1006-2 was de- 
rived by subtracting the Lya flux from the total flux in the 
NB1006 magnitude. This Lya deficit might be caused by the 
gradual and continuous reionizing process from z = 5.7 to 7.3, 
because the density of the rest-UV continuum is insensitive to 
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FIG. 10. — Number density jilvq, and star formation rate density SFRDl vq of the LAE at z = 7.215, SXDF-NB 1006-2 confirmed spectroscopically in this study 
(filled square) and those at 2.3 < z < 7 from the spectroscopic studies down to L Lya = 1 .0 X 10 43 erg s . The filled circles are those of the LAE at z = 6.96 (Ota 
et al. 2008). The crosses show the densities at z = 5.7 and 6.6 in the SDF calculated from the samples of (Kashikawa et al. 2011). Densities at 2.3 < z < 4.5 
and z ~ 4.5 are calculated using the best-fit Lya Schechter LFs from van Breukelen et al. (2005) (inverse triangles) and Dawson et al. (2007) (pentagons). Also, 
densities at z = 3.1,3.7,5.7 (diamonds) and 6.6 (triangles) in the ~ 1.0 deg- of SXDF are calculated using the best-fit Lya Schechter LFs from Ouchi et al. 
(2008, 2010). Each horizontal error bar shows the redshift range of each survey. The vertical error bars at z = 5.7,6.6,7 and 7.3 include both cosmic variance 
and Poissonian errors for small number statistics. We correct for the detection completeness of the z = 7.3 LAE survey. If SXDF-NB1006-1 is a real LAE at 
z = 7.3 with a brighter L Lya than one of SXDF-NB 1006-2, the number and SFR densities of SXDF-NB1006-2 becomes small the open square slightly shifted in 
the redshift direction for clarification. 
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FIG. 1 1 . — Cumulative rest-UV continuum LFs at z = 5.7 (blue curve) and 
6.6 (red curve) derived the LAE samples in the SDF (Kashikawa et al. 201 1). 
The filled circle and square is the rest-UV number density of the z = 6.96 (Ota 
et al. 2008) LAE and the SXDF-NB 1006-2, respectively. The UV magnitude 
of the SXDF-NB 1006-2 is comparable with that of the rest-UV bright end at 
Z = 5.7 and 6.6. 



the surrounding HI clouds and might not evolve. Note how- 
ever that the UV magnitude of SXDF-NB 1006-2 has a large 
error caused by Lya flux error (see Table 4), and the Lya 
deficit might arise from the galaxy evolution as well as HI 
clouds. This tendency is consistent with the results of the pre- 
vious surveys for LAEs z < 7 (Ota et al. 2008; Kashikawa 
et al. 2006, 2011; Ouchi et al. 2010) and our results might 



show for the first time that the trend continues up to z ~ 7.3. 
On the contrary, our result is inconsistent with the results of 
Hu et al. (2010); Tilvi et al. (2010); Krug et al. (2012); Hi- 
bon et al. (2010) which reported that Lya LFs do not evolve 
from z = 5.7 to 7.7. If the decline of Lya luminosity is caused 
only by IGM attenuation, the neutral fraction xhi of hydro- 
gen gas at z = 7.3 can be estimated to be ~ 0.79 using a ra- 
diative transfer simulation based on McQuinn et al. (2007) 
which calculate the number density of LAE at the epoch of 
reionization, and we assume that the IGM at z = 5.7 is com- 
pletely ionized. Compared with L Lya of SXDF-NB 1006-2 and 
the Lya LF at z = 7.3 constructed by the model of Kobayashi 
et al. (2010), the IGM transmission for Lya photons T[™ can 

be estimated. Here, T[™ is calculated with T[™ = observed 

L Lya I model L Lya . As we use the faintest L Lya in the low- 
est number density within their errors as the observed L Lya , 
the derived T[™ is given as an lower limit, 0.27. Thus, the 

possible T[™ ranges 0.27 < T[™ < 1. If we apply the cal- 
culated T[™ to the Lya flux attenuation model of (Santos 
2004), the neutral fraction xhi at z = 7.3 is estimated to be 
0.82. Note that the derived neutral fraction is strongly model 
dependent. However, it must be noticed that one remaining 
candidate, the SXDF-NB 1006-1, has yet not been observed 
spectroscopically. If this candidate was a real LAE with sim- 
ilar Lya luminosity with the SXDF-NB1006-2 at = 7.3, the 
number density would increase to that at z = 7, and not evolve 
between z = 7 and 7.3 as indicated by small points in Figure 9 
and 10. On the other hand, if all the candidates were variable 
objects, which indicates that no LAE at z = 7.3 were detected 
in our survey volume, the derived xhi is only given as an upper 
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limit. 



6. SUMMARY AND CONCLUSIONS 

We have carried out a systematic imaging and spectro- 
scopic survey for LAEs at z ~ 7.3, for the first time, in the 
SDF and the SXDF, using the Suprime-Cam equipped with 
new red-sensitive CCDs and a newly developed NB1006 fil- 
ter. The total survey area and volume are 1719 arcmin 2 and 
5.94 x 10 5 Mpc 3 , respectively. Four z ~ 7.3 LAE candidates 
were photometrically identified from the NB1006 imaging 
down to the NB limiting magnitude of 24.6-24.8, and we ob- 
tained spectroscopy for three of them, SDF-NB 1006-1, 2 and 
SXDF-NB 1006-2 using Subaru/FOCAS and Keck/DEIMOS. 
We have detected a definitively asymmetric emission line 
from SXDF-NB 1006-2. If this line is Lya, it is a LAE at 
2 = 7.215. SXDF-NB 1006-2 has a luminosity of 1. 2^ % x 10 43 
erg s -1 , a star formation rate 11 M Q yr" 1 estimated from the 
Lya emission line, and a weighted skewness S w =4.90 ±0.86. 

The conclusions of this survey are summarized below. 



• Compared with the number density and SFR density 
of LAEs at z < 7, The number density and SFR den- 
sity of LAEs monotonically decreases from z = 5.7 to 
7.3 while they do not change from z = 3 to 6, and 
the number density of rest-UV continuum at z = 7.3 is 
comparable with the rest-UV cumulative LFs at z = 5.7 
and 6.6. The derived decreases of the number den- 
sity are n% 7 a 3 /n z = 5 a 7 = 0.05!^ from z = 5.7 to 7.3 and 

«i;"Kya 6 = °- 16 ±ai6 from z = 66 t0 13 down t0 
j_Lya = J Ox 10 43 erg s" 1 . This trend is consistent with 
the results of the previous surveys for LAEs at z < 7, 
and it might be caused by IGM attenuation. 

• The magnitude of SDF-NB 1006-2 in the NB1006 im- 
ages varies and it might be a AGN/QSO at z = 7.3 or 
1 .7 rather than a LAE, in spite of the detected asymmet- 
ric emission line (S w = 4.59 ± 1.25). The third object, 
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show any significant variation, though the time span to 
check the variability in SXDF was not long enough to 
completely trace the typical variability. Hence, it is im- 
portant to the check whether the candidates are transient 
objects or not by dividing photometric images into sev- 
eral epochs. 

In order to accurately measure the neutral fraction and re- 
veal the relationship of LAEs and LBGs in the evolution of 
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generation facilities such as the Subaru/Hyper Suprime-Cam 
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